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ABSTRACT 

Wc construct barotropic stellar equilibria, containing magnetic fields with both poloidal 
and toroidal field components. We extend earlier results by exploring the effect of dif- 
ferent magnetic field and current distributions. Our results suggest that the boundary 
treatment plays a major role in whether the poloidal or toroidal field component is 
globally dominant. Using time evolutions we provide the first stability test for mixed 
poloidal-toroidal fields in barotropic stars, finding that all these fields suffer insta- 
bilities due to one of the field components: these are localised around the pole for 
toroidal-dominated equilibria and in the closed-field line region for poloidal-dominated 
equilibria. Rotation provides only partial stabilisation. There appears to be very lim- 
ited scope for the existence of stable magnetic fields in barotropic stars. We discuss 
what additional physics from real stars may allow for stable fields. 

Key w^ords: 



1 INTRODUCTION 

At least three different classes of star appear to harbour strong and broadly similar magnetic fields: Ap/Bp stars on the 
main sequence, magnetic white dwarfs and neutron stars. All three classes exhibit long-lived fields, with large-scale simple 
structure. Surface fields are up to the order of 10* G for Ap/Bp stars, 10^ G for magnetic white dwarfs and lO'^^ G for 
neutron stars. Although these values are very different, the resultant total magnetic fiux is similar for the three classes, due 
to their correspondingly different radii. Given their high strengths , magnetic fields are expected to p l ay important roles in the 



evolution and dynamics of these stars (IDonati fc Landstreet 



200£ 



addition, magnetic fields distort a star ( Chand rasekh ar fc Fermil igBSl) : this distortion mav (in a rotating neutron star) lead 



Wickramasinghe fc Ferrariolboool : iHarding fc Lajbood ). In 



to appreciable emission of gravitational waves (Bona zzola fc Gourgoulo n 1996). 

Observations only give us information about the external field of a star, with no probe of the interior. This is a problem, 
since the details of the field geometry can greatly affect the ratio of internal-to-external field strength; the observed field may 
provide a poor estimate of the global value. For example, if a star's magnetic field is confined to an outer region (like the 
crust of a neutron star), its volume-averaged value could be lower than expected from the observed surface value; if the field 
is predominantly buried within the star, its strength may be higher than expected. In numerable stu dies have been motivated 
by the attempt to understand the kinds of magnetic field that may exist in stars (see Mestel ( 19991 ) for a review). 

One approach to the modelling of stellar magnetic fields is to construct equilibrium configurations of a fiuid star, and 
then to show that they are dynamically stable (and hence credible models for the long-lived magnetic fields present in many 
real stars). In practice, the latter step is very difficult — it requires confirmation that every possible perturbation about 
a stellar equilibrium is stable. For the former step, constructing hydromagnetic equilibria, virtually all studies make the 
simplifying assumption of a barotropic (i.e. unstratified) stellar model: the pressure is taken as a function of densit y alone. 
This has been criticised as unrealistic — temperature/entropy profiles provide stratification in main-sequence stars ( MesteJ 
19561 ) and white dwarfs jReiseneggeilboogl). whilst the varyin g proportions of charged particles (composition gradients) result 
in neutron-star stratification ( Reisenegger fc Goldreichll 19921 ). Stratification provides an additional degree of stabilisation, and 
numerical simulations for main-sequence stars with entropy gradients have provided strong evidence for the existence of stable 
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magnetic equilibria in this case jSraithwaite fc Nordlundlbood : iBraithwaitj I2OO9I I . It is not clear if this result is applicable to 
neutron stars, however, where the origin of the stratification is different. 

Despite its limitations, a barotropic equation of state is a sensible first approximation to stellar matter, and one which 
is astrophysically relevant if it allows for stable magnetic equilibria. It is still unknown whether such stable configurations 
exist, but there seem two reasons to expect them to. The first is connected with the geometry of a mixed poloidal-toroidal 
field. A poloidal field is known to be unstable in the region where its field lines close within the star, whilst a toroidal field 
suffers localised instabilities around the pole; in a mixed field, each field component 'fills in' the region where the other 
field compon ent is unstab l e. This chan ges the magnetic field geometry in these unstable regions, potentially suppressing the 



instabilities (|Wrightlll973 



for barotropic equilibria ( Lander fc JonesI 20091: 



Tavlei 198C). The second reason to ex pect stable equilibria is that the 'twisted-torus' fields found 



Ciolfi et al. 



2OO9I) are qual itatively similar to the final state of the nonlinear 



stratified-star simulations mentioned above ( Braithwaite fc Nordlund 200a). 

This paper aims to explore the possibility of stable magnetic fields existing in barotropic stars, by testing the stability of 
a number of equilibrium models. We begin, in section 2, by exploring the range of magnetic equilibria possible in barotropic 
stars, reviewing the key equations of the problem and examining the various classes of solution. Broadly, the resulting 
equilibria differ in how their magnetic field and current are distributed. In section 3, we discuss properties of perturbations 
in a magnetic star: the governing equations, how a mixed poloidal-toroidal field breaks the equatorial symmetry of the 
system, and potential instabilities. Next, we use the mixed-field equilibria of section 2 as background configurations for time 
evolutions of the linearised perturbation equations (section 4). We study perturbations known to result in instabilities for 
purely poloidal/toroidal fields, to determine whether or not these perturbations are still unstable in mixed-field stars. Following 
this, we discuss some important aspects of stellar physics not accounted for in barotropic models, and whether these may 
remove magnetic instabilities (section 5). Finally, in section 6 we summarise our results. 



2 MIXED-FIELD BAROTROPIC EQUILIBRIA 

The first step towards studying the stability of magnetised stars is to construct suitable ba ckground equilibrium models, which 
we describe in this section. Although full details may be found in iLander fc JonesI liooi), we pause here to focus on the key 
equations of the problem, and any restrictions imposed by our assumptions. This is important, since our aim is to produce 
the widest possible range of equilibrium models to test for stability. 

We model a star as a barotropic fluid body, in Newtonian gravity and axisymmetry. The star has a mixed poloidal-toroidal 
magnetic fleld B, and our scheme also allows for rigid rotation (about the magnetic symmetry axis) with angular velocity 
Q.. We work in cylindrical polar coordinates (117,0,2), with the z-axis being aligned with the symmetry axis of the star. An 
equilibrium will be described by the stationary form of the MHD Euler equation: 

where P is fluid pressure, p mass density, $ gravitational potential and j the electric current. This equation needs to be solved 
together with Ampere's law 

47rj = V X B, (2) 



Poisson's equation 

and the solenoidal constraint on the magnetic fleld 



4^Gp, (3) 



V ■ B = 0. (4) 



The system is closed with an equation of state. We choose a polytrope with index N = 1, & typical approximation for neutron 
star matter: 

P = P{p) = = kp\ (5) 

Note that a more suitable value for white dwarfs would he N = 1.5, whilst main seq uence stars tend to be modelled as A'^ = 3 
polytropes. We are able to generate magnetic equilibria for a variety of values of N ( Lander fc Jone j 20091 ) . but the resulting 



field configurations are all qualitatively similar. Since magnetic-field instabilities are related to the geometry of the fleld (see 
section [373)) . we have some expectation that our results for = 1 should be representative of the stability of main-sequence 
stars and white dwarfs, not just neutron stars. 

Some algebra is needed to recast the equilibrium equations above into a convenient form for solution. First let us take 
the curl of the Euler equation (U]), which — in the case of a barotropic star — yields 

= VM, (6) 
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for some scalar function M. The solenoidal constraint on the magnetic field together with the assumption of axisymmetry 
allows us to express the poloidal field component Upoi in terms of a streamfunction u: 

Bpoi = — Vu X e^. (7) 



Furthermore, the toroidal component is related to the streamfunction through some function /: 

B4> = -f{u) 



(8) 



and it may also be shown that A'l = M{u). So, the magnetic field and the force it exerts on the fluid are both related to the 
single scalar function u. Equivalently, it is possible to work with the 0-component of the magnetic vector potential A instead, 
since u = wA^. With vari ous algebraic tricks, one may derive a relation between these functions u,M and /, known as the 
Grad-Shafranov equation ( Grad fc Rubin 19581 : Shafranov 19581 ): 



Au 



f 



du 



1 d 



+ ■ 



w dvj dz 



(9) 



Combining the Grad-Shafranov equation with Ampere's law ([2} in axisymmetry yields a useful relation between the magnetic 
functions, the current and the fleld: 

AM 



J 



47r Au 



- pvj- 



Au 



(10) 



By rewriting j and B in terms of the streamfunction ijTomimura fc Eriguchill2005l l we arrive at a version of Poisson's equation 
for the magnetic field: 



V ( — It sin ) — 



f Af , AM\ 
zu Au Au J 



(11) 



The above result is the last that one can obtain without loss of generality. For numerical solution we now need an integral form 
of equation (|lip . and in choosing a Green's function to produce this integral we also implicitly choose a boundary condition 
for the exterior field. All of the equilibria we generate use the integral equation 



-u sm — 



1 

4-K 



f'{u)f + 'iTvmpM'{u) 



' df , 



(12) 



for which the magnetic-field magnitude falls off as a dipole, B — > 1/r'^ as r ^ oo. If one wished instead to solve for fields 
confined within the star (for example), a different Green's function would be required. 

One must also specify the functional forms of M{u) and /(u), which are chosen on physical grounds. In particular, we 
see from equation (|10p that the derivatives M'{u) and f'{u) are related to the distribution of electric current. Both M and / 
can be chosen to allow for smooth or discontinuous current distributions inside the star, whilst most choices for / result in 
currents outside the star. As discussed later, if one wishes to avoid exterior currents the choices for /(u) become very limited. 
The integral equations for $ and u are solved together with the Euler eq uation using a numerica l schem e which iteratively 



finds equilibrium configurations — for details see iLander fc Jones! (|2009l l or lTomimura fc Eriguchil (|2005l l. We now consider 



different choices of these functions and the resultant equilibria. 



2.1 Exterior poloidal field, no surface currents; 'type 1' 

Our first choice of surface treatment is that employed in iLander fc Jonesl (|2009l ). Specifically, the toroidal component is 
confined within the star and goes to zero smoothly at its surface, whilst the internal poloidal field matches smoothly to an 
external component, which falls off as We feel this is the most natural boundary condition for a fiuid star in vacuum, 

although it ignores the different physics present in the outer regions of stars — for example the crust, ocean and magnetosphere 
of a neutron star. One surprising, perhaps unsatisfactory, aspect of these solutions is that we find the toroidal component is 
always weak with respect to the poloidal one in a global sense; its contribution ftor to the magnetic energy fmag is never more 
than a few per cent of the total. Locally the two field components can be more equal, however, with the maximum values of 
poloidal and toroidal components being of comparable magnitude. 

To produce the desired surface behaviour of the toroidal field component, we need to make a suitable choice of magnetic 
function f{u). The only functional form which achieves this appears to be 

f^i^U ^tsurfmax) ^ ^ ^surfmax 
U ^ ^surfmax 

is the maximum surface value attained by u (jTomimura fc Eriguchil I2OO5I : ILander fc Jonesl I2OO9I : ICiolfi et al 



(13) 



where Usur 



200^. This choice has the effect of enclosing the toroidal field component within the closed-field line region of the star, 
producing equilibria where the toroidal field only occupies a small volume of the star. Nonetheless, it is the largest volume 
that does not give rise to an exterior current. 

We consider this choice slightly inconsistent with the requirement that / be a function of it alone: Usurfmax is a constant, 
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type la type lb typc Ic 

0.15 I , , , , , 1 0.12 I , , , , , 1 0.06 I , , , 




-0.15 I 1 1 1 1 1 1 -0.08 I— ^ — - - - - I - --- ---r -•" 1 1 1 1 -0.1 I 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 

Figure 1. Three different type 1 or 'twisted-torus' equilibria, all with a relatively large toroidal component: £tor/^'mag = 3%, 4.5%, 4.5% 
for types la, lb and Ic respectively. The variation in magnitude of the three field components is shown as a function of dimensionless 
stellar radius r/Rt, where Rt is the surface of the star. Bg and are evaluated along the equatorial axis, whilst Br (zero along this 
axis) is evaluated along the polar axis. The neutral line of the poloidal component is located where the Bg line crosses the x-axis, and 
the toroidal component is centred around this line. These three equilibria represent the range of solutions we have found which have an 
exterior field but no surface/exterior current; nonetheless, all are qualitatively similar. 



but one found by evaluating m at p = 0. The value of ttsurfmax varies between iterative steps, and hence the scheme involves an 
implicit density dependence. Instead, it is possible to produce entirely consistent equilibria using the streamfunction's global 
maximum itgiobmax (which is attained within the star): one can fit the / function to contours of different fractions of Wgiobmax 
and pick the solution with the largest volume of toroidal component which does not extend outside the star. The resulting 
equilibria are virtually identical to those using Usurfmax, however, suggesting that (|13p is in fact an acceptable choice. 

All of our type- 1 equilibria use the above choice for /, but differ in the choice of the other magnetic function M. In 
Lander fc Jones (|2009l ) we found that only M = ku {k = const) produced magnetised equilibria, with other cases resulting 



in our numerical scheme iterating to unmagnetised solutions. Here we refer to the M = ku equilibria as type la. 

More recently, we have found a way of using other functional forms of M without the numerical scheme iterating to a 
zero-field solution. By including a 1/Wgiobmax factor, the field at each iterative step is 'normalised' to the same magnitude as 
the previous step and not driven to zero. These different functional forms of M produce qualitatively similar equilibria, so we 
take one representative example — M = fcii^/ugiobmax — and refer to it as type lb. 

A different generalisation of the original type-la equilibria is to relax the requirement that the current be continuous 
inside the star. Allowing for steps in the current means we can choose forms of M and / with discontinuous tt-derivatives. We 
find that a step in f'{u) produces little qualitative difference in equilibria from types la or lb, so we do not consider these 
further. More significant differences emerge by choosing a step in M' (u) — as an example of these we choose the following: 

M'(^) = |'' ">'^-'t— (14) 

10 ti ^surfmax ■ 

Note that this corresponds to a step in the interior current, but not to a surface current; see equation H10[). Since we wish to 
avoid 5-distribution behaviour of M'{u) (and hence the Lorentz force) in this class of solutions, however, we must ensure that 
M itself is smooth — and so we take 

X «:(" — Wsurfmax) U>U 

M{u) = { (15) 

I U 5^ ^surfmax- 

Note that this means the Lorentz force is zero outside the closed-field line region, but that a force-free field can still exist in 
the rest of the stellar interior. The other magnetic function / is chosen in the same way as for type la. This final choice we 
refer to as type Ic. 

The three type-1 field configurations are compared in figure [T] where we plot the variation in magnitude of Br,Bg and 
with stellar radius. The maximum value attained by the toroidal field is greater in type lb than type la, but it is still 

confined to the same small region. The profiles of types la and Ic differ mainly in the relatively weak radial field in type Ic. 

Overall, however, the three field configurations appear very similar. 



2.2 Confined fields, no exterior poloidal component; 'type 2' 

Next, we consider a very different class of equilibrium from that of the previous subsection: one where the magnetic field is 
confined to the stellar interior. We will refer to such confined-field equilibria as 'type 2'. As described earlier, the boundary 
condition on the magnetic field is imposed by the choice of Green's function used to solve equation to generate equilibria 
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0.15 




r/R^ 

Figure 2. A type-2 equil i brium , where all field components are confined within the star, not just B^. Plotted here is the lowest-eigenvalue 
solution of lHaskell et al.l bOOSi ). corresponding to no nodes in the toroidal field. 65% of the magnetic energy is contained in the toroidal 
component — far higher than for type-1 equilibria. 



with confined fields we would need to change this function. Since the total magnetic field vanishes at the stellar surface in 
this case, there would be no separate need to ensure that the toroidal field vanishes outside the star. 

The major disadvantage in assuming a confined field is the more limited astrophysical relevance: in particular, many stars 
have strong exterior magnetic fields, thought to be comparable in strength with the internal value. One possible application 
of this treatment could be to model accreting neutron stars or white dwarfs, where the surface field is weak but a strong 
internal field could still exist, 'buried' under the accreted matter. 

There are interesting differences between the properties of magnetic equilibria in the confined and non-confined cases. 
The non-confined models described above (type 1) are 'poloidal-dominated' — the poloidal field component contains most 
of the star's mag netic en ergy — whilst confin ed-field solutions seem to be toroidal-dominated by the same measure (jloka 
200 ll : iHaskell et al..,200a : iDuez fc Mathislboiol ') . Related to this, type-1 magnetic fields induce oblate distortions in the star's 



density distribution, whereas type-2 equilibria have prolate ellipticities. Finally, since the magnetic fields in the two cases are 
qualitatively different, it is natural to expect them to have different stability properties. 

Generating confined-field equilibria would not be a straightforward extension to our work, nor is there strong motivation 
to do so from observations. On the other hand, we would like to compare t he two contra s ting fi eld geometries of type 1 and 2 
equilibria. To this end, we turn to the perturbative semi-analytic models of lHaskell et al. I (|2008l ) for confined mixed-field stars. 
These authors found that the confined-field case reduced to an eigenvalue problem, with higher eigenvalues corresponding to 
a larger toroidal component with an increasing number of nodes. We take their lowest-eigenvalue solution for the magnetic 
field together with a spherical density distribution, thus neglecting the (small) distorting effect of the magnetic field. Unlike 
our type-1 equilibria, this is not fully self-consistent, but the resulting models should be accurate enough for a qualitative 
idea of stability. 

We plot the magnetic field we are adopting as 'type 2' in figure (2] Even for this lowest-eigenvalue solution, the toroidal- 
component magnetic energy is 65% of the totaB Clearly, the toroidal field is far more globally significant than for type-1 
equilibria, extending throughout the stellar interior. 



2.3 Exterior poloidal component, surface currents; 'type 3' 

The types of magnetic field outlined in the last subsections represent two extremes of mixed-field equilibria with a confined 
toroidal component; in type 1 the poloidal component is non-zero across the surface and into the exterior, whilst in type 2 
it vanishes at the surface. For both types, the toroidal component goes to zero smoothly at the surface. Type 1 equilibria 
are always dominated by the poloidal component; type 2 by the toroidal component. Since we wish to test the dependence 
of stability on the relative strength of the two field components, we would like a way of introducing a 'sliding scale' between 
types 1 and 2, to produce equilibria with more equal proportions of each field component. In addition, the type-1 boundary 
treatment neglects the different nature of the outer region of a star, where the mass density becomes very low and we can expect 
differences from ideal MHD, due to the increasing resistivity and (for a neutron star) the crust, ocean and magnetosphere. 

As a crude attempt to investigate the effect of different physics in the outer regions, and as a way to produce mixed- 
field equilibria with stronger toroidal components, we consider equilibria with surface currents in this subsection. We stress 
that we are not asserting that surface currents themselves are necessarily significant in real stars — they simply provide a 

^ Tabic 1 from iHaskell et al ] l|2008l l indicates that Stov/Smag = 90%, due to a normalisation error: their averaged poloidal field is 
~ 2.2 X 10^2 G, not 1 x lO^^ G as reported. 
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r/R* 

Figure 3. Equilibria with the toroidal field component matched to a surface current rather than going smoothly to zero there. For the 
type 3a model, £tor/£'mag = 50%; for type 3b it is 30%. The right-hand plot is an attempt to produce an equilibrium with a surface 
current that also looks roughly like a twisted-torus structure. 



mathematically convenient way of producing a surface boundary condition which is different from that of the type 1 or 2 
equilibria discussed above. 

For these equilibria we choose a form of / such that is non-zero up to the stellar surface. An exterior toroidal field 
requires an exterior current; since we still want to avoid this we appeal to some (poloidal) surface current to match the non-zero 
interior toroidal field to a zero exterior one. This surface current is not separately modelled; its form is whatever required to 
produce a co nsistent matching for t he toroidal field component. A similar approach to the construction of magnetic equilibria 



was taken by Colaiuda et al. (200y). In practice, then, only the form of the magnetic function / is changed; 



f , -. \ Cb interior 

/'H = r , . (16) 

(J exterior 
I au interior 

/(«) = , . (17) 

I exterior. 

We refer to this choice as type 3a. Note that since f(u) — zuB^, this corresponds to a jump in the toroidal field component 
at the stellar surface. 

The toroidal field is now distributed throughout the stellar interior. As a half-way choice between this type of equilibrium 
and the usual twisted-torus structures of type la (where the toroidal field only exists in some small torus near the surface), 
we consider one final class of equilibrium — type 3b — designed to look like a twisted-torus configuration but with a larger 
torus than in type 1. Again, we use a surface current to achieve this: 

j'(^) _ I" interior and it > 0.5tisurfmax ^^g^ 
I otherwise 

I a(u - 0.5usurfmax) interior and u > O.Susurfmax 

= i „ , . (19) 

10 otherwise. 

The relative magnitudes of the r, 9 and 4> components of the magnetic field in these surface-current equilibria are plotted 
in figure [3] Type-3a equilibria look quite similar to the confined-field solution in figure O differing mostly near the surface, 
since the surface- current equilibria have exterior fields. The type-3b equilibrium was constructed to resemble a twisted-torus 
(i.e. type-1) equilibrium, but with a larger toroidal component; cf. figure [T] 

We conclude this section by comparing the various magnetic equilibrium models we have constructed. In figure |4] we plot 
poloidal-field lines and represent the toroidal component magnitude with coloured shading. Since the three type-1 equilibria 
are rather similar, we include only one of these (type la) as a representative example. 



3 PERTURBATIONS OF A MIXED-FIELD STAR 
3.1 Perturbation equations and numerics 



Given the apparently intractable nature of the analytic problem ( Tavlei 1980l l . we will study the stability of mixed poloidal- 
toroidal field stars using a numerical approach. Our code evolves the linearised MHD perturbation equations, taking the 
equilibria of the previous section as background configur ations. We have successfully used this code to stud y instabilities in 
stars with purely toroidal fields ( Lander fc Jonej 2011a ) and purely poloidal fields ( Lander fc Jonej 2011bl ). and so expect 
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type la type 2 




Figure 4. Comparison of various magnetic-field equilibrium structures. The poloidal field lines are shown in black, whilst the strength 
of the toroidal field is represented by the colour code. The black arc from a; = 1 to y = 1 is the stellar surface. Type la is a typical 
'twisted-torus' configuration, where the toroidal field component sits inside the closed field lines. The other field configurations do not 
have this twisted-torus geometry. Type 2 is a confined-field, where all poloidal-field lines are closed, whilst types 3a and 3b are equilibria 
with surface currents. 



its results for mixed fields to be reliable. Since full details of our numerical scheme are reported in lLander et al. 

I l|2O10l ). we 

content ourselves with a brief summary of the salient details here. 

We wish to study the behaviour of linear perturbations about MHD equilibrium; this allows us to test the stability of the 
equilibrium, based on whether or not we find unstable modes of the system. We work in the Cowling approximation, where 
perturb ations in the gravitat ional potential are not evolved. This approximation is known to overestimate the stability of a 
system jMoss fc Tavleill 19691 ). placing some potential doubt over any magnetic fields we find to be stable, but not over those 
systems we find to be unstable. 

Instead of working directly with the velocity v and the perturbed magnetic field SH, we define fiux variables f = pov and 
(3 = poSJi (zero-subscripts denote background quantities). This allows for a simpler, more numerically stable treatment of the 
stellar surface. Our final perturbation variable is the perturbed density Sp. Now, in the frame corotating with the background 
star, the linear perturbations of the system are described by the following time-evolution equation^: 



— = — V(5pH — 2f2xf 

at po po 

--^(V X Bo) xBo + — !— (V X Bo) XI3+ -^(V x /3) x Bq - -^(Vpo x /3) x Bo, (20) 
47rpo 47rpo 47rpo 47rpo 

^ = Vx (f xBo)-^ X (f xBo). (22) 
ot Po 

This 3D system of equations may be reduced to 2D, using the axisymmetry of the background to decompose in the 
azimuthal angle <j). We then run evolutions for a chosen azimuthal index m, which is convenient for studying the stability of 
particular modes. Imposing the perturbations' symmetries at the poles then reduces our numerical grid to one half-disc. For 



^ Note that in our earlier papers using this time evolution code, algebra errors affected the second and third terms of the perturbed 
Euler equation. The equations used in the code were, however, correct. 
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a purely poloidal or purely toroidal field, equatorial symmetries allow a further reduction of the grid to a single quadrant of 
a disc, but these symmetries are broken for a mixed poloidal-toroidal field; see the next subsection. 

The code evolves equations using a predictor-corrector algorithm and is second-order convergent. We include artificial 
viscosity — a fourth-order Kreiss-Oliger dissipation term added to the Euler equation — but take care to do so at the 
minimum value required to damp numerical instabilities, so physical instabilities are minimally affected. No corresponding 
artificial resistivity term was added to the in duction equation. T o prevent growth of V ■ B errors, we employ an auxiliary 
variable and equation for divergence-cleaning (jPedner et al.ll2002l ). 



3.2 The importance of equatorial symmetry 

The equatorial symmetry of perturbations is closely related to the behaviour of modes and instabilities in magnetic stars. 
Here we explore this important aspect of the perturbation problem and introduce some useful notation for later. 

Let us start with the behaviour of perturbations of an unmagnetised fiuid star. For this case, the four perturbation variables 
split into two sets based on their shared symmetry properties; we denote these sets 5i = {fr, fif,, Sp} and 5*2 = {fe}- For 
example, an initial perturbation in /,. which is symmetric (or 'even') about the equator will induce corresponding symmetric 
perturbations in and 5p and an anti- symmetric fe perturbation. A different initial perturbation, leading to equatorially 
antisymmetric behaviour of the perturbations in class 5*1, would produce symmetric perturbations in the 5*2 element. 

From this, we can define two symmetry classes for modes (coherent global responses of the fluid). The first we refer to 
as , and corresponds to 5*1 elements being even and S2 elements odd; the other class, , is that of odd Si elements and 
even 5*2 elements. Note that the class is equivalent to the set of modes variously called polar, polar-led or spheroidal — 
for example, the /-mode. The class is the set of modes known as axial/axial-led/toroidal, the most familiar example being 
the r-mode. 

If we now add a purely poloidal magnetic field to the background, an analysis of the equations shows that the original 
81,82 classes are augmented by magnetic variables in the following manner: 5*1 = {fr,f(p,Sp,l3g} and 5*2 = {/s, /3r, /3</)}. The 
chief instability for a purely poloidal field is a kink mode in the closed field line region, as described in the next subsection. This 
corresponds to motion in the 8 direction across the equator and hence to a symmetric perturbation there; an antisymmetric 
ve is zero along the equator. For this reason the instability occurs only in the "P" symmetry class, and its development can 
be prevented in evolutions which enforce "P^ conditions at the equator. 

In the case of a star with a purely toroidal magnetic field, the division of perturbations based on equatorial symmetry 
changes from the poloidal-field case and the two sets become 81 = {fr,f4,,Sp,l3r,P4,} and 5*2 = {/s,/?e}. Since Si and S2 
take different forms in the purely poloidal and purely toroidal cases, it follows that when the background field is mixed 
poloidal-toroidal, the equatorial symmetry of the perturbations will be lost. As a result, we will no longer have a division into 
■p^ and classes, and instead the two will be 'coupled' through the magnetic field. Equally, there will no longer be distinct 
axial-led or polar-led modes in the conventional use of the terms. 



3.3 Diagnosing instability 

Both purely poloidal and purely toroidal magnetic fields are susceptible to rapidly-developing instabilities, of essentially the 
same origin: the unstable nature of a cylindrical plasma column. Temporarily neglecting gravity, the simplest instabilities of 
the cylindrical system are the sausage or varicose mode and the kink mode, illustrated in figure [5] (left-hand side). Based on 
this, one can identify regions of a star where magnetic instabilities may occur: the cylindrical region around the pole in a 
toroidal field, and the torus around a poloidal field's neutral line (the line along which the poloidal field vanishes). 

We now turn to the right-hand side of figure [5] and consider first a star with a purely toroidal magnetic field. The 
geometry is similar to that of the cylindrical plasma, but the star's self-gravity must now be accounted for. Since both sausage 
and kink instabilities in this case can operate parallel to shells of constant gravitationalpotential, however, neither should 
be suppressed; accordingly, we expect m = and m = 1 magnetic instabilities to exisl|j. In a poloidal field the geometry 
is different, with the unstable cylinder closed into a torus, resulting in potential instability in all m > perturbations. The 
sausage mode of the cylinder, on the other hand, is suppressed by gravity; so we do not expect a significant m = instability 
for a poloidal-field star. 

For a mixed poloidal-toroidal field to be stable, we need to check that all instabilities that occur for a purely poloidal or 
toroidal field are eliminated in the mixed configuration. In this paper, we only consider non-axisymmetric modes. To test for 
any instabilities due to the toroidal component, therefore, we look at the behaviour of m = 1 perturbations in the vicinity of 
the pole. By contrast, we anticipate that instabilities of the poloidal component will be localised in the closed-field line region. 



^ In the special case of an incompressible star, the sausage mode does have to move fluid elements along the z-axis and hence will be 
suppressed by gravity, but the m = 1 instability will remain. 
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Figure 5. Left: instabilities of a cylindrical plasma column, with field lines in the 0c-direction labelled with arrows. The two main 
instabilities, in the absence of gravity, are the sausage/ varicose mode (with cylindrical azimuthal index rric = 0) and the kink mode 
(rric = 1). Right: the corresponding potentially unstable regions of a mixed-field star — the cylindrical region around the pole where the 
toroidal component vanishes, and the torus around the neutral line, where the poloidal component vanishes. The potentially unstable 
values of m, now accounting for gravity and working in the spherical polar coordinates of the star, are shown. 



for any m ^ 1. For a purely poloidal field only initial data will produce these (kink) instabilities, but a mixed field allows 
for p+~'P~ coupling — so that any initial data may result in the excitation of a poloidal-component instability. 

In the initial phase of an instability, an unstable mode is excited with an amplitude which grows exponentially in time; 
to study this, it is sufficient to consider the behaviour of linear perturbations about an equilibrium. Beyond some amplitude 
the mode 'saturates', causing a nonlinear rearrangement of the magnetic field for which the perturbative regime is no longer 
applicable. Since we only evolve the linearised system, however, the background is stationary and experiences no reduction in 
magnetic energy as the amplitude of a perturbation grows — hence the instability is able to grow indefinitely. This is clearly 
not realistic, but it does allow us to determine the growth rate of the initial instability with ease. In all cases, the onset time 
of the instability is expected to correspond to the time taken for an Alfven wave to cross the relevant part of the system. As 
a diagnostic, we define the following quantity as an order of magnitude estimate for this time: 

_ i?. 

where ca is the Alfven speed and overbars denote volume- averaged quantities. We take a field strength of B = 3 x 10^'' G as 
our canonical value, for which the dimensionless Alfven timescale ta ~ 40. This strong magnetic field and correspondingly 
short dynamical timescale allows for faster numerical evolutions, but we have confirmed our results also hold at lower field 
strengths. In addition, we have checked that the growth rate of any unstable mode converges with numerical grid resolution, 
and scales approximately linearly with the magnetic field strength. 

Note that confirming the stability of a magnetic field is far more difficult than showing it is unstable: for the latter 
one need only find a particular mode with exponential growth in time, whilst for the former one must establish that every 
potential perturbation results only in innocuous oscillations about the equilibrium and not unstable growth. Within the 
limitations of our numerics, we are able to test the stability of equilibria against nonaxisymmetric initial perturbations with 

or V~ symmetry, u p to m — 6. Since this test easily shows up instabilities of purely poloidal and purely toroidal fields 



y or y symmetry, u p to m — b. bmce tins test easily snows up mstabuities or purely poloidal and purely toroidal nelds 
( Lander fc Jonejboilal lbh. we may have some confidence that a mixed poloidal-toroidal field showing no such instability is 



'stable'. More accurately, it will not be susceptible to the fastest-growing class of instabilities and is likely to be dynamically 
stable; we cannot establish its stability on secular timescales. 



4 STABILITY ANALYSIS OF MIXED-FIELD STARS 
4.1 Type 1 (twisted-torus) equilibria 

We now begin our analysis of the magnetic stability of barotropic stars containing both poloidal and toroidal components, 
using the various background models discussed in section 2. In particular, we are testing for the possible presence of any 
instabilities known from the purely-poloidal and purely-toroidal field cases. In this subsection we test the stability of twisted- 
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poloidal-ficld background mixed-field background 




i 

Figure 6. The instability of type-1 (twisted-torus) equilibria; results are shown for type la, but the corresponding plots for types lb 
and Ic are very similar. We compare the behaviour of m = 1 and 2 initial perturbations (with and V~ symmetry) in a star with 
a purely poloidal field (left), and in one where a toroidal component has been added (right). For the mixed-field star ftor/£^mag = 3%; 
a high percentage for this class of equilibrium, but not producing any stabilisation. For a mixed field, magnetic perturbations across 
the equator have no definite equatorial symmetry; in a sense the mixed field allows for 'mixing' of the two 'P-classes. For this reason, 
although initial data starts by evolving stably, as for a poloidal-field background, at later times energy is transferred to the unstable 
class of perturbations. 



torus equilibria (type 1), as constructed in section [2.11 These equilibria are poloidal-dominated, with a toroidal component 
only occupying the small volume within the closed-field line region. 

In figure|6]we compare the stability of a type-la mixed field with that of a poloidal field. It would b e reasonable to expect 
the addition of a toroidal component to remove, or at least reduce, the instability of the poloidal field ( Wright 19731 ): despite 
this, our results show no such stabilisation. On the left-hand side we plot the evolution of perturbed magnetic energy for 
P"*" (polar) and 'P~ (axial) perturbations, on a poloidal-field background. As expected from section the class evolve 
stably, whilst the V class suffer an instability that appears at t ~ 30 — 50; cf. our Alfven timescale estimate of ta ~ 40 for 
this field strength (B = 3 x 10^*^ G). 

The right-hand plot shows the corresponding evolutions for a mixed-field star with ftor/i^mag ~ 3%, relatively high for 
a type-1 equilibrium. As for the poloidal field, after one Alfven timescale the perturbations grow exponentially; furthermore, 
since the mixed-field breaks the equatorial symmetry of the perturbations (see section |3^ . even initial data results in an 
excitation of unstable modes. The growth rate of the m = 2 instability is not reduced by the addition of a toroidal component, 
nor are the m = 4 or 6 instabilities (omitted in figure|5]to avoid cluttering it). The m = 1 instability, in fact, grows faster for 
a mixed field than a purely poloidal one. Finally, although we have only plotted results for type-la equilibria, type-lb and Ic 
mixed fields are similarly unstable. 

To confirm that the instability from figure [6] is due to the poloidal component, we look at its location within the star. In 
figure [7] we plot the magnitude of the perturbed field |(5B|, in the case of a poloidal field and a type-1 mixed configuration. 
The plots are in logarithmic scale and show the late stages of the evolution, where the unstable mode is completely dominant. 
For both field configurations shown, there is clearly exponential growth around the neutral line — slightly more localised for 
the mixed field. Although we only show results for m = 2, the unstable growth occurs in the same location for m = 1. In this 
latter case we see no additional growth around the pole, which would have been indicative of a toroidal-field instability. If 
such an instability does exist in this case, its growth rate must be much slower than the dominant poloidal-field instability. 



4.2 Type 2 (confined-field) equilibria 

Next we turn to a stability analysis for a star with a confined magnetic field, as described in section [2^ This model does not 
have a twisted-torus structure and is not poloidal-dominated; 65% of its magnetic energy comes from the toroidal component. 
These differences allow us to explore whether the twisted-torus structure or dominant poloidal field are responsible for the 
instability of type-1 equilibria. 

On the left-hand side of figure [8] we plot the evolution of perturbed magnetic energy for m — 1,2,4 and 6 azimuthal 
indices. Although the m = 2, 4 and 6 instabilities still appear to exist for this magnetic field, their growth rates are greatly 
reduced in comparison with the type-1 equilibria results. In addition, the onset time for the instabilities is rather longer: about 
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Figure 7. The magnitude of the perturbed magnetic field |(SB|, in logarithmic scale, at late times in the evolution of two stars. The 
left-hand plot is for 'P~ perturbations on a poloidal-field background, whilst the right-hand plot shows perturbations of a mixed field 
(for which there is no equatorial symmetry); m = 2 for both cases. The poloidal-field instability — manifested as exponential growth 
around the neutral line — is seen to operate in the mixed-field star too. The corresponding plot for m = 1 shows growth around the 
neutral line too, but none around the pole, so we have no evidence for additional instabilities originating from the toroidal component 
of the mixed field. 
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Figure 8. Left: the instability of a confined-field (type 2) equilibrium, seen from the exponential growth of the perturbed magnetic 
energy in time. The background stellar model is very different from type 1; it is toroidal-dominated (ftor/^^mag = 65%) and does not 
have a twisted-torus geometry. Although the field strength is the same as before (B = 3 X 10^^ G), no instabilities manifest themselves 
until t ^ 200 — five times our estimated Alfvcn timescale. Right: the unstable growth in |(5B| is fastest around the pole, suggesting that 
the toroidal component is the origin of the instability for type-2 equilibria. 



five times that for type 1. One similarity between the type-1 and type-2 results is that the m — 1 instability seems to grow 
more quickly for a mixed field than a purely poloidal one. 

Since the type-2 equilibrium is toroidal-dominated, we wish to see which field component is responsible for its instability. 
Accordingly, on the right-hand side of figure [8] we plot |5B| throughout the star for an m = 1 evolution. We see the fastest 
unstable growth occurring near the polar surface, indicating that the toroidal component is now unstable (see section l33|l . 
There is also substantial growth away from the pole, however, and m > 1 perturbations are seen to be unstable too; both of 
these results indicate that the poloidal component also plays a role in the instability of this field type. 
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4.3 Type 3 (surface-current) equilibria 

In this subsection we test the stability of type 3 equilibria, our final class of mixed-field configurations, where we have an 
exterior poloidal field but nonetheless a strong toroidal component. This is achieved by allowing for a step in the toroidal 
component at the surface and requires a corresponding surface current. Recall, however, that we are not able to study the 
surface dynamics of the perturbed magnetic field, since we evolve the fiux variable f3 — poSJi. 

We consider two examples of surface-current equilibria, type 3a and type 3b, described in section[2]3] Both have substantial 
toroidal components, with comparable energy to the poloidal component: for type 3a ftor/fmag = 50%, and for type 3b 
ftor/£^mag = 30%. In cach case, the value was the upper limit on the percentage of toroidal field possible with our chosen 
functional forms. From figure |9] we see some degree of stabilisation in these equilibria, compared with those of type 1. In the 
evolution for type 3b, the m — 6 instability appears to have been removed and the m = 4 instability considerably reduced. 
However, the m = 1 and 2 instabilities remain. 

Type 3a is the closest to a stable barotropic equilibrium we have found, with no evidence for the existence of m = 2, 4 
or 6 instabilities. Despite this, we still see exponential growth for m — 1 — which is enough to render this magnetic field 
unstable. Given the stabilisation of other modes, one might question whether the m — 1 growth is numerical rather than 
physical in origin. We have checked this, as for all other unstable configurations, by convergence-testing the instability growth 
rate. Furthermore, the unstable growth was localised around the neutral line, closely resembling the plots in figure [T] This 
suggests an instability of physical origin, related to the poloidal component of the magnetic field, even though the toroidal 
component contributes 50% of the magnetic energy. 



4.4 Rotation 



We conclude our stability analysis by looking at the effect of rotation. Rotation may act to reduce magnetic instabilities, 
though there is some uncertainty in the literature about how effective the stabilisation will be. Part of the problem is that 
there are no conclusive ana l ytic results for the rotating case, due to the increased comple xity of the stability conditions 
Frieman &: Rotenberj 1960l : Lvnden-Bell fc Ostriker 1967 ; Glampedakis fc Andersson 2007 ) with respect to the static case 



Bernstein et al. 



19581 ) 



Although no definitive analytic results exist for the stability of rotating magnetised stars, the work of lPitts fc Tavleij (|l985l ) 
provides some suggestions, based on approximate treatments of some model problems. They predict that rapid rotation may 
remove magnetic instabilities already present in the non-rotating star, but that a new instability will be introduced, albeit 



field instabilities ( Braithwaite 20061 : Kitchatinov fc Riidieei 



necessarily provide complete stabilisation. For p oloidal fields . 
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that rotation has some stabilising effect, whilst IBraithwaitd (|2007l ) does not. With somewhat tentative expectations, we now 
turn to our results for a rotating star with a mixed poloidal-toroidal field. We specialise to m = 2 perturbations on a type-la 
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Figure 10. Left: the evolution of tlie perturbed magnetic energy in a rotating mixed-field star, for different dimensionless rotation rates. 
Rotation is seen to reduce, but not remove, the instability. Note that the Kcplcrian frequency in these units Q,k ~ 0.72. Numerical 
instability prevents us from studying more rapidly-rotating stellar models than Q = 0.412, around 57% of the Keplerian value. Right: 
instability growth rate as a function of rotation rate, for two field strengths. The instability present for f2 = seems to be reduced by 
rotation, but at around Q = 0.15 the lines plateau, suggesting the appearance of a second instability. 



mixed-field background in this subsection, but the results are representative of other field types and perturbations we have 
tested. 

To study the effect of rotation on magnetic instabilities, we monitor the evolution of the perturbed magnetic energy once 
again. We quantify the (exponential) growth rate of instability using the parameter 
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(24) 



From the left-hand plot of figure [10] we find that rotation certainly decreases the growth rate of the perturbations, although 
some degree of unstable growth is present in all results. Our code suffers some numerical instabilities for very rapidly-rotating 
stars, so the high est rotation rate we i nclude here is 57% of the Keplerian value. 

The work of lPitts fc Tavleil ijlQSSh suggests that a non-rotating star s magnetic instabilities should be suppressed when 
the rotational angular velocity exceeds the corresponding Alfven velocity Qa ~ 2tt/ta- Using our previous estimate for ta 
(f23l). we expect this suppression to occu r at » » 0.11 for B ^ 2.25 x lO^'' G and at 0.14 for B = 3 x 10^® G. As the 
original instability is removed, however, iPitts fc Tavler (|l985l ) predict a new instability will be induced whose growth rate 
may be slower but still significant. 

With these results in mind, we now turn to the right-hand plot of figure 1101 Here we quantify the effect of rotation, 
plotting the rate of exponential growth of the magnetic energy as a function of rotation rate. We consider stars with B — 2.25 
and 3 x 10^® G (a ratio of 3:4), to investigate the scaling of instability growth rate with field strength. At zero rotation rate 
the more highly magnetised star has a growth rate 4/3 that of the other, confirming the linear dependence we expect in 
this case. Adding rotation causes an initial decrease in growth rate which is roughly linear in $1, up to some threshold value: 
fl « 0.13 for the weaker magnetic field, Q « 0.16 for the stronger field. These values agree rather well with the estimates of 
the previous paragraph. Beyond the threshold rotation rate for each star, the unstable growth rate seems to reach a plateau 
value of ~ 10% that of the non-rotating case. We believe this corresponds to the new magneto-inertial instability predicted 
to occur for sufficiently rapid rotation. 

5 DISCUSSION: POSSIBLE STABLE EQUILIBRIA 

The aim of this paper has been to construct as wide a range of mixed poloidal-toroidal magnetic fields as possible in barotropic 
stars, and th en to test their stability. We have found that all of these equihbria are unstable: poloidal-dominated twisted- 
torus mo dels ( Tomimura fc Eriguchil 2005 ]_ Lander fc JonejboogI : Iciolfi et al.ll2009l ) , toroidal-dominated fields confined within 
the star (|Haskell et allboosl : lOuez fc Mathijboiol ) and equilibria with surface currents jColaiuda et ahlboOijl ). For equilibria 
where ftor/fmag ^5 50% the instability's origin appeared to be the poloidal component; for the confined-field model we tested, 
with £tor/fmag ~ 65%, wc fouud cvideucc of a toroidal-field instability, but some indications of instability in the poloidal 
component too. This seems to leave little scope for configurations where all magnetic instabilities are suppressed. On the 
other hand, we have observational evidence for long-lived stellar magnetic fields, so stable magnetic equilibria clearly exist. 
Here we discuss possible mechanisms to suppress the instabilities we have found. 
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Rotation: highly magnetised stars (like many Ap/Bp stars, magnetic white dwarfs and magnetars) often rotate relatively 
slowly. In some, however, the kinetic rotational energy greatly exceeds the magnetic energy and one might expect suppression 
of magnetic instabilities. We find, however, that when rotation seems rapid enough to suppress the original magnetic instability 
(of a non-rotating star), a new instability appears to be introduced. This new magneto-inertial instability has a growth rate 
of ~ 10% the non-rotating value — still very short on stellar timescales. 

Non-axisymmetry: our code uses a decomposition of the perturbation equations in the azimuthal angle, which relies 
on having an axisymmetric background star. It is possible that a n on-axisymmetric m agne tic field will not s uffer from the 
same instabilities we have found. Indeed, the non-linear evolutions of lLaskv et al.l pOllh and lCiolfi et al.l (|201ll ') seem to show 
unstable initial fields developing non-axisymmetric structure. At the end of these reported simulations, however, the field 
remains highly dynamic in the closed-field line region, so it is probably too early to regard the results as stable equilibria. It 
would be very interesting to see whether these fields settle down over longer evolution times. 

Stratification: perhaps the easiest way to explain our results is to say that they are irrelevant for real stars, where 
the equation of state is not barotropic but instead stratified. Indeed, the only numerical evolutions which seem to produce 
stable magnetic equilibria use a stratified stellar model, where the stratification is due to temperature/entropy gradients 
jBraithwaite fc Nordluiidl bood : iBraithwait j boogj ). Whether these equilibria are, in fact, stable is not completely settled; 



Bonanno fc UrpinI (|201ll ) discuss a generic magnetic instability whose lengthscale can be very short, and hence not show up 



in numerical simulations (for which the grid is necessarily relatively coarse). 

In any case, results for entropy-based stratification may not be applicable to neutron stars, where t he stratification comes 



from composition gradients. Recently, some first attempts have been made to model such equilibria ijMastrano et al.l 12011 



Glampedakis et al.l l2012l : iLander et al.l l2012l ) , but these are qualitatively similar to barotropic results and their stability is 



unknown. Note that there is no guarantee that ra dial stratification will suppress th e instabilities we find, which mainly involve 
6-direction motion (although it does appear to in iBraithwaite fc Nordluiidl (|2006l ') and related papers). 

Elastic crust: in the case of neutron stars, the outermost layer of the star (around a kilometre in thickness) consists 
of an elastic crust. This crust does have the potential to suppress unstable motion of the kind studied in this paper, unless 
the magne tic field is so strong that its energy exceeds the crustal elastic energy; this occurs when the field strength exceeds 
~ 10" G ijThompson fc Duncanlll993h . although the precise value depends on the poorly-known yield strain of the crust. 
Beyond the critical value, a magnetic instability can cause the crust to crack, which is a plausible scenario for the triggering 
of magnetar bursts ([Thompson fc Duncanlll995l ). For white dwarfs and Ap/Bp stars, there is no analogous outer region to 
hold an unstable field in place. 

Superconductivity: relatively little is known about the equilibria and stability of magnetic fields in superconducting 
stars. The magnetic force takes a very different form from the norm al Lorentz force, however, so it is not u nreasonable to 
expect significant differences in the magnetic field configuration s too ( Mend«l| 19981: Glampedakis et al. 20 111 ) . Neutron stars 



quickly become cool enough to contain superconducting regions ijPage et al.ll201ll : 
physics may help to stabilise (or destabilise!) their magnetic fields. 



Shternin et al.ll201ll ) , so this aspect of their 



6 SUMMARY 

We construct a wide range of barotropic stellar equilibria, with mixed poloidal-toroidal magnetic fields. These include poloidal- 
dominated and toroidal-dominated configurations, as well as those where both components have comparable energies. Every 
field configuration is found to suffer instabilities whose origin appears to be one of the two field components. Starting with 
a poloidal-dominated field, increasing the strength of the toroidal component suppresses unstable poloidal-component modes 
of higher azimuthal number, but the m = 1 instability seems difficult to remove. For a sufficiently strong toroidal component, 
the origin of the m — 1 instability changes from the poloidal to the toroidal component. Rotation reduces the instability 
growth rate, but seems to introduce a new instability of slightly weaker growth rate (around 10% of the non-rotating value). 

Although this paper does not prove the non-existence of stable fields in barotropic stars, it suggests that any such fields 
form — at best — a very restricted class (which we have been unable to obtain). By contrast, many classes of star are observed 
to have long-lived magnetic fields, so we expect that stable magnetic equilibria do exist. Our results indicate that barotropic 
models are, therefore, of limited astrophysical relevance. For main-sequence stellar models, the addition of entropy-gradient 
stratification seems to allow for stable equilibria, and the same may then be true for white dwarfs. For the modelling of neutron 
star fields, we believe that it is now particularly important to account for additional physics beyond the barotropic-fiuid model, 
including the effect of a crust, composition gradients and superconductivity. 
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